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Abstract 


This paper introduces the inverse dissipation potential into the Newtonian Dynamic 
equation and studies the motion equations of the objects in the isolated gravitational 
system. It is found that at large scales it can derive the dynamical equation of cosmic 
expansion similar to the ACDM model and yield the flat rotation curves for spiral 
galaxy. Different from the usual dark matter models, the derived flat rotation curves 
are the result of time accumulation rather than the direct action of mechanics. And the 
Tully-Fisher relationship is also discussed, it is found that the basic constant ao in the 
MOND model and the form of the function u have a clear corresponding physical 
significance in the model of this paper. 
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1.Introduction 

Dark matter and dark energy are two major problems in today's cosmology [1-4]. 
These two "invisible" substances show two total distinct effects that dark matter 
produces gravity, while dark energy produces anti-gravity. Various theoretical models, 
based on general relativity, are proposed to explain these two phenomena. Obviously, 
there are two different modification methods. By changing the left side of the Einstein 
field equation, the time-space geometry itself is modified and form the theories 
represented by /(R) gravity theory [5,6], brane-world gravity theory [7], MOND [8,9], 
etc. The other is to modify the right side of the Einstein field equation with additional 
fields or matter within the time-space, such as ACDM [10], Quintessence [11], 
phantom [12,13], etc. 

However, none of the present theories could solve all problems perfectly, 
especially, since the search for dark matter particles and the measurement of dark 
energy have found nothing [14,15], it prompted us to think that perhaps modifying the 
gravity equation was an option to try. And it can be seen in Ref. [16-21] that some 
scholars have done a lot of work in this area. 

Among the various models that replace dark matter, the MOND model proposed 
by M. Milgrom is a popularly discussed one. This model presented an empirical 
formula, in which a basic constant ao is introduced to explain some observational 
phenomena [22-28], for example, it can explain the Tully-Fisher relationship well 
[24,25]. However, the model still has some problems in explaining other observation 
data, such as the mass deviation in the cluster of galaxies [29], and the theory itself as 
it violates the basic law of momentum conservation [30]. 

In short, the MOND theory is not perfect at the time of its submission. Therefore, 
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based on the original MOND theory, Bekenstein et al., TG Zlosnik et al. and M. 
Milgrom proposed the relativistic MOND theory, namely Tensor-Vector-Scalar 
Theory (TeVeS) [20,26], Einstein-Aether theory [31] and subsequently Bimetric 
MOND Theory [32,33], respectively. Although these theories greatly enrich the 
content of MOND, problems still exist [34,35]. 

The partial success of MOND theory and its extension theories vaguely implies 
that it maybe stems from a more basic theoretical system. Hence, this paper, based on 
the theory of non-relativistic Newtonian dynamics, utilizes the theoretical method of 
correcting inertia [37] by introducing the dissipative potential to Newton's Dynamic 
equation to discuss its influence on the motion of objects at large scales, and attempts 
to further reveal the deep physical significance of the basic constants ao in MOND 
theory. 


2.Newtonian dynamic equation with dissipative potential energy 
In the frame of non-relativistic, the Newtonian dynamic equation is 
F=ma (1) 
Where m is the mass of the object and a is the kinematic acceleration. 

There are usually two ways to modify the Newtonian dynamic equation: the first 
is to modify the right side of the equation to establish a modified inertia theory, and 
the second is to modify the left side of the equation to establish a modified gravity 
theory. In the case of relativity, the two modified theories are equivalent. While, in 
non-relativistic case, and applying the modified inertia theory, Eq. (1) can be assumed 
to be the following modification form 

F=m(a-Ay) (2) 
Where J is a constant, v is the velocity of a moving object. 

For an isolated gravitational system 2, consisting of only two objects, and an 
object with mass m is moving within a 2D-plane in the gravitational field of a 
quasi-stationary object with mass M (M>>m), we can obtain the Lagrangian equation 
for system » based on Eq. (2) 


r=- = m roje 


(3) 


Where (n 0) is the polar coordinate system of the 2D-plane, G is the 
gravitational constant. 

As a velocity-dependent dissipation function is introduced in Eq. (2), then the 
corresponding dissipation potential is 


P= -im +(r0)] (4) 


According to the Hamiltonian principle of a dissipative system, the motion 
equation of the object with mass m is 
Le 0)=2r° 0 
dt 
eo 2 GM . (5) 
r-r0 +—=Ar 
r 
When the object with mass m moves along the radial direction of the polar 
coordinate system, namely 0=0. Thus, the second formula of Eq. (5) can be rewritten 
as 
r=- ~ +Ar (6) 


2 
F 


As we known that the cosmological equations can be obtained from Newtonian 
mechanics [37]. For a sphere with the radius of r, if the physical coordinate r=a(t)R 
(R is the comoving coordinate of the cosmic expansion, a(t) is the cosmic expansion 
factor), the density is p=pva%, and the quality of the material inside the sphere is 
M=4z/3pr’, it can be expressed from Eq. (6) that 


ee 4 e 
ey, (7) 
3a 
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Substitute Hubble's law H = al a into Eq. (7), we can obtain 
_ 4aGpy 
3a’ 
Comparing Eq. (8) to the time-time component of the field equation in the 
ACDM model [10], it can be obtained 


+ AHa (8) 


A 
| AH (9) 
Eq. (9) shows that the Newtonian dynamic equation with inverse dissipation 
potential can derive the dynamic equations similar to the expansive cosmological 
ACDM model at large scales. 
Now we discuss the second case. For a given initial condition at t=0, that an 


object with mass m circle a quasi-stationary object with mass M in its gravitational 
2 


field, ie 1 =0, GM/r, =0, 1. Ina dissipative system, r and 0 will change with 


time. Let’s assume r= fit), 0 = 0, f(t), then the second formula in Eq. (5) can 
be rewritten as 


.. e2 1 GM . 

ITA tar = Ar (10) 

“ GM 1 . 

r+-——(—-fifr)=Ar (11) 
rm f 


Considering the initial condition, when the gravitational field is very weak, 
namely 
GM 


Bai (12) 
k 
Then substituting Eq. (12) into Eq. (11) 
PaA (13) 
Solving Eq. (13), then the radius r is 
r=re° (14) 
From the first formula in Eq. (5), it can be obtained that 
rO=r 6, e” (15) 
Therefore 
rô =r, (16) 


As the rotational tangential velocity V =6r , the initial rotational tangential 
velocity V, = 6, n , then 
V =V; (17) 


It can be indicated from Eq. (14) and Eq. (17) that the radius increases 
exponentially with time, but the rotational tangential velocity of the object remains 
the same, which coincides with the asymptotical flat rotational velocity properties 
typically exhibited by some spiral galaxies [38,39]. 


3.MOND Model 

As we konwn, M. Milgrom introduced a basic constant ao into the MOND theory, 
which modified the actual acceleration and Newtonian acceleration in Newtonian 
mechanics as follows [8,9] 

ay =au(a/ ay) (18) 

When a/ao>>1, H (a/ao)~1; when a/ao<<1, H (a/ao)~a/ao. 

From Eq. (12), it shows that when the gravitational field of the initial condition is 
very weak, the tangential velocity of the object with mass m remains basically 
unchanged with time, while, the distance r could still increase. Suppose the intensity 


of the weak gravitational field at initial corresponds to go, that is 
GM 


&=— (19) 
k 
Substituting Eq. (19) into Eq. (17), it can be obtained 
V*=Vi = GMg, (20) 
Rewriting Eq. (20) 
lg(L)=4lg(V)-lg(Gg,<M/L>) (21) 


Where M/L is the mass-to-light ratio. 

We can see that Eq. (21) is just the Tully-Fisher relationship. Comparing Eq. (20) 
with the MOND theory, it can be found that go=ao. Therefore, 
the centripetal/actual acceleration of the object obtained from Eq. (14) is 

y? Vy Eo _ 4 
a (22) 


r ne e e 


The Newtonian acceleration is 
_GM_ GM By . Gy 


ay 2.2 24 2M 2A (23) 
F rje e e 
Divide Eq. (23) by Eq. (22) 
a 1 a 
2A = Tht — (24) 
a e ay 


Hence, comparing Eq. (24) with Eq. (18), it shows that the model in this paper 
can derive the form of the empirical function u introduced in the MOND theory. 


4.Conclusion 

An isolated gravitational system is studied in this paper through introducing the 
inverse dissipation potential into Newton's dynamic equation. It is found that the 
modified dynamic equation can simultaneously deduce the cosmic expansion dynamic 
equation similar to the ACDM model and the asymptotical flat rotational velocity 
properties typically exhibited by spiral galaxies, which is generally considered to be 
dark energy and dark matter, respectively. It suggests that the dark energy effect and 
the dark matter effect may be two aspects of the same physical reality after the 
introduction of the inverse dissipation potential. 

Unlike the usual dark matter model, this kind of the asymptotic flat rotational 
velocity properties of the spiral galaxies are the result of long-term accumulations 
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with time, rather than the direct action effect of mechanics. And in its derived cosmic 
expansion dynamic equation, an equivalent time-varying cosmological constant is 
presented, which is consistent with the Quintessence model and the phantom model. 

Finally, same with the MOND model, the modified Newton's dynamic equation 
also derives the Tully-Fisher relationship. Moreover, the basic constant ao introduced 
in the MOND model and the form of the function u show a clearly corresponding 
physical significance in the model of this paper. However, the discussion in this paper 
is still not a relativistic correction of the Newtonian gravity, which requires further 
study. 
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